Potential of a moving test charge in a dusty plasma in the 
presence of grain size distribution and grain charging dynamics 



Michael A. Raadu*and Muhammad Shafiq^ 
Royal Institute of Technology, Alfven Laboratory, 
Division of Plasma Physics, SE-100 44 Stockholm, Sweden 

February 2, 2008 

Abstract 

It is well known that the form of grain size distribution strongly influences the lin- 
ear dielectric response of a dusty plasma. In previous results [IEEE Trans. Plasma 
Sci. 29, 182 (2001)], it was shown that for a class of size distributions, there is 
an equivalence to a Lorentzian distribution of mono-sized particles. The electro- 
static response to a slowly moving test charge, using a second order approximation 
can then be found [Phys. Lett. A 305, 79 (2002)]. It is also well known [Phys. 
Plasmas 10, 3484 (2003)] that the dynamical charging of grains in a dusty plasma 
enhances the shielding of a test charge. It seems natural at this stage to seek the 
combined effects of grain size distribution and grain charging dynamics to a test 
charge moving through the dusty plasma. Here we consider the effects of both grain 
size distribution and dynamical grain charging to a test charge moving slowly in 
a dusty plasma by expressing the plasma dielectric response as a function of both 
grain size distribution and grain charging dynamics. Both analytical as well as the 
numerical results are presented. It is interesting to note that the previous results 
can be retrieved by choosing appropriate values for different parameters. This kind 
of study is relevant for both laboratory and space plasmas. 
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1 Introduction 



There is currently considerable interest in understanding the physics of dusty plasmas 
which in addition to the electrons and ions, also contain a dust component. The dust 
grains become charged due to the collection of ions and electrons from the plasma and can 
typically acquire thousands of electron charges (due to high mobility of electrons, dust 
grains usually becomes negatively charged). Moreover, the dust charge does not necessar- 
ily remain constant and may continuously fluctuate due to varying plasma currents that 
flow onto the dust charge surface. The currents reaching the dust grain surface depend 
on the ambient plasma conditions and the floating potential of the dust particle. In this 
way the dust charge becomes a dynamic variable and fluctuates about an equilibrium 
charge state. In order to deal with the problem of charging dynamics, many papers have 
taken into account this dynamics and presented their results [U 121 IH] , the consequences 
of including a dust component have lead to a renewed interest in the problem of the test 
charge response. This is important for understanding the influence of a dust component 
on the interaction between charged particles. An important consequence of the potential 
excited by a moving charge is the energy loss and braking of the velocity due to the 
resultant electric field at the moving charge [U E] . 

We have investigated the response of a slowly moving test charge in a dusty plasma 
in the presence of charging dynamics (HI [Zj and found that the dynamical charging of 
dust grains in a dusty plasma enhances the shielding of the test charge. The response 
potential up to the second order in test charge velocity was found and expressed analyti- 
cally in terms of strength functions. A delayed shielding effect due to dynamical charging 
was also reported. The linearised dielectric theory was used and the equilibrium dust 
distribution was considered to be Maxwellian. Furthermore, the equilibrium dust parti- 
cles were assumed to be similar and all the dust particles were treated as point particles. 
But this is not always true and in general, a size distribution of dust grains is to be 
expected both in artificial and natural plasmas [HJ |§] , for a Maxwellian distribution and 
a special class of physically reasonable size distributions, the dielectric response function 
was shown to be equivalent to that for monosized particles with a generalised Lorentzian 
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or kappa distribution (THj Recently, we have taken into account the test charge response 
of a dusty plasma with a grain size distribution |TT) IT2] and have shown that the form 
of grain size distribution strongly influences the linear dielectric response of a test charge 
in a dusty plasma. The analytical expressions for the response potential, using a second 
order approximation were found and the effects of collisions also investigated. More re- 
cently, A. M. Mirza et al. ^3] extended this work further and presented analytical as 
well as numerical results for the slowing down of a pair of test charge projectiles moving 
through a multicomponent dust-contaminated plasma. In their analyses, they found that 
the energy loss for the Maxwellian distribution is larger compared to that for generalised 
Lorentzian distribution. They also found that for smaller values of the spectral index k, 
the test charge projectile gains instead of losing energy. 

In this paper, we have extended our previous work on grain size distribution by 
taking into account the effect of charging dynamics and present analytical as well as 
numerical results for the response potential for a test charge moving through a multicom- 
ponent dusty plasma. 

2 Plasma Dielectric for a Dusty Plasma 

The linear response of the dusty plasma for an electrostatic disturbance can be determined 
through the choice of the plasma dielectric function. Here the dielectric will include a 
term for the dynamical charging of the dust grains and the effect of a specific choice for 
the size distribution will be taken into account. 

2.1 Grain Size Distribution 

Here we choose the size distribution h(a) used previously [Hi]. 



where the constant ho is defined by setting the integrated density to the dust density 
n d (h = 3n d a l3+1 /V{{(3 + l)/3)). The distribution h(a) has a maximum at a = a = 
(P /S) 1 / 3 ^ 1 . If we let a — ► oo with (3 ~ a 3 the distribution h(a) tends to a delta function 




(1) 
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at a = do i.e. a monosized distribution with dust grain radius do- This limit is useful for 
comparing the general results that will be found here with earlier results for monosized 
distributions. The distribution h(a) can also be transformed to a distribution over grain 
mass m so that, 

h(a)da = w(m)dm = wom^~ 2 ^ 3 exp (— /xm) dm (2) 

For small sizes h(a) has an approximate power law dependence on the size a. A power law 
dependence is a simple first approximation if the actual size distribution is not known. 
For large masses w(m) is dominated by an exponential decrease with mass (for (3 = 2 the 
dependence is purely exponential, as assumed in early work on interstellar dust grains). 
These properties motivated this particular choice of size distribution ^Uj. Without charg- 
ing dynamics this choice for h(a) leads to a dielectric response equivalent to a kappa 
distribution [TU] with k = (2/3 + 5)/6. (For the purely exponential mass dependence, 
/3 = 2, the index k = 3/2.) 

2.2 Charging Dynamics 

Here we now include charging dynamics with frequencies u = Q u0 and Q v0 (defined by 
Melands0 et al [T^]) that depend on grain size a. This leads to a response term with an 
integration over grain size with a differential "charging wavenumber" H^ch given by 

Hdch ( a ) ^ a = 47r— ^-a/i (a) da 

where Q u q and Q v o are the frequencies introduced by Melands0 et al ^H] in the linearized 
equation for the grain charge perturbation g^i with a plasma potential perturbation 1; 

dq d i 



dt 



-tt u0 q dl - 47Te aQ V Q(j)i (3) 



In this equation Q u q = vq{o) is the grain charge relaxation rate i.e. tq = 1/uq is the time 
scale for the grain charge to come into equilibrium with the undisturbed plasma. The 
total charging wave number K dch may then be defined by 

oo 

K 2 dch = I Hl h da 
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For a monosize distribution h(a) is a delta function and the above expression reduces to 
the standard definition [7j. Integrating for h(a) given by the equation (0) gives for H dch 
the expression 

K dch = 47i— — -A— (4) 



fi u0 a r (2±l) 



2.3 Plasma Dielectric 



For a general size distribution with charge relaxation rate vo(a) that is a function of the 
dust grain radius a the plasma dielectric function is, 



K 2 K 2 



K 2 
' K 2 



2k \ (K-V t \ fK-Vt 



2k-V\ V td J V« 



1 ?_ff|Ja)^o(a) 



o 



where .fTo arid are the effective Debye wave-number and effective thermal velocity for 
the dust as defined in jTHj. For V t < Vtd, the plasma dispersion function Z K [K ■ V t /V t dj 
is given as follows [TB] : 

/ at-va = ^ - MV" r(. + f(n + 2) ) ^-VA" f6 , 

"V V td J KV 2 T(K-\)tA W r(i(n + 2)) \ V td ) K> 
Following the analysis of Melands0 et ai ^H] , for a standard model of the dust charging 
process, explicit expressions can be found for the frequencies Q u0 = u (a) and Q v0 . These 
may be written as [7j, 

&vo — fivo\ — u P i (7) 

*De 

&u0 — ^uOT ^pi (8) 

*De 

where, assuming equal ion and electron temperatures, the numerical constants are 5 v q = 
2.793 and 5 u q = 1.795. Here we note that, for grain sizes comparable to the electron 
Debye length X De , these frequencies are of the order of the ion plasma frequency u pi . 
The frequencies Q v o and f2 u0 = are simply proportional to the dust size a and the 
ratio 5o = Q v o/fl u o = 8 v q/8 u q = 1.556 is independent of the dust size. The last term in 
equation © can now be written using these expressions and the size distribution h(a) 
defined by equation ffl- There is no obvious simple analytical expression for the resulting 



integration, but for a slowly moving test charge the integral can be expanded as a power 
series in V t . The individual terms can then be integrated in terms of the gamma function. 



3 Response to a Moving Test Charge 

For a test charge response in a plasma, the general expression for the electrostatic potential 
is given by [Hj 

Qt f exp[.K.r] (g) 



r 8ir 3 e J K 2 D(K } K • V t j 
where Vt is the test charge velocity and D(K,u) is the plasma dispersion function. The 
explicit form of D(K,u) depends on the physics of the dusty plasma. Here D(K,uj) is 
chosen to include the effects of a grain size distribution and charging dynamics. 

For a slowly moving test charge (Vt < V t d), we can expand the plasma dispersion 
function (equation ©) up to first order and hence rewrite equation (0) for the dielectric 
up to second order in test charge velocity as 



K 2 D(K,K-V t ) K 2 + K 2 



eff 
2 



1 + iA (0) 



K 2 D fk-v t 



K 2 + K 2 ff \ V td 



K 2 + K* ff \ V td J o u0 uj pi K 2 + K'ff v 

°u0 u pi + n eff 



-1 



(10) 



with the definitions 



r ( l3+2 ) 

K\ = 4rr5 ^, K* ff = K 2 De + K 2 Di + K\ + K 2 f ((3) ,/(/?) = 



and with 



M + i T (1/3 + i)' 2/3 + 5 ' r(ttl) 



3^ ^ 6 V3^ 1 3J V 3 y 

Here the term K 2 f(j3) appearing in the definition of K e ff is identical to Kj ch given by 
equation (j3J). The total charging wave number Kd c h measures the contribution of charging 
dynamics to the total effective shielding wave number K e ff. In the above relations f3 is the 
power law index of the size distribution for small radii, related to the equivalent kappa 
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distribution by the relation k = (2/3 + 5)/6. Expanding equation (fTUJ) for the inverse 
dielectric function up to second order in test charge velocity V t and using in equation (1), 
we may express the electrostatic potential as <fi = 0i + 4> c h with 



Qt 



exp[iK ■ r] 
7T% J K 2 + K 2 e 



B{(3)K 2 D fk-V t 



eff 
2 



1 - 



iA{(3)K 2 D (K-V t 



K2 + K tff 



K 2 + K 2 eff 



V 



Id 



(K 2 + K, 



2 

eff 



V v td 
k-v t 



dK 



and 



Ich 



Qt 



exp [iK ■ r] 



+C{(3) 



7T% J K 2 + K eff L 
« 2 AL K\K 2 
5 2 u0 u 2 pi K 2 + K 2 f 



—i 



a\ De K 2 K ,~ 
5 u0 uj pi K 2 + K 2 ff y ' * 
2 a 2 X 2 De K\K 2 



-2A 09) 



aA De K K 2 D K 2 



(k.v t ) 



X2 , ,2 
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'2 



K-V t 



(K 2 + 



2 



(#-v t ) : 



(12) 



<5«0 U pi V td (]£2 _|_ ^2 

It is to be noted that is the same as we found earlier ^2] except for the definition 
of K e ff which now includes the effect from charging dynamics in terms of K\, while <p c h 
is the contribution which comes explicitly from the dust charging dynamics. The reader 
is referred to J2| for the results of equation (fTTj) for 1; while in the following we shall 
present the results for cf)^. The above equation (|T2*|) can be written in terms of strength 
functions as 



(f>ch(r, A) 



qt 



87T%r 



V t ffu(r) cos A + V t 2 (g 20 (r) + g 22 {r) cos 2 A) + 0(V t 3 



(13) 



where A 
strength 

9ii 0) 
920 (r) 



is the angle between the test particle velocity V t and the radial vector r. The 
functions 9ij{r) are given by the following expressions 



aA 



7T 



De 



8u0 ^pi 



K 2 exp(-rK eff ) 



4K eff ^0 U 2 pi 



*K 2 ff C (P) 
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- 1 
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eff 



exp {—rK, 



eff) 



* vtd °uO U p i r J\ e ff 

- (r 2 K 2 eff + 3 + rK eff ) (r 2 K 2 ff + 3 - rK eff ) * + 2rK eff (9 + r 2 K 2 eff ] 

where the following relations defining $(rK e ff) and $?(rK e ff) in terms of exponential 
integrals [7j have been introduced, 

$ (rK eff ) = exp (rK eff ) E 1 (rK eff ) - exp (-rK eff ) E { (rK eff ) (14) 

* (rK eff ) = exp (rK e// ) ^ (rK eff ) + exp (-r^ e// ) ^ (rK eff ) (15) 

$(?/) and \l/(|/) (for y > 0) are directly defined as principal parts of integrals (here, for 
real y, equivalent to taking the real part) as follows, 

ny)=Re[-r 2te f-f dt ) (16) 







2 exp(— yt) dt 



* j < 17 > 

From these definitions it follows that $(?/) = d^(y)/dy and that ^(y) = d&(y)/dy + 2/y. 
The functions $(?/) and \I/(?/) introduced here are closely related to the auxiliary functions 
f(y) and g(y) used in the analysis of the Sine and Cosine Integrals [T7j. As for f(y) and 
g(y) asymptotic forms may be found for Q(y) and ^(y), 

- i * (y) ~ 2/~ 2 + 3! y- 4 + 5! y- 6 + 7! y" 8 + O (y" 10 ) (18) 

| * (y) ~ y- 1 + 2! y~ 3 + 4! y" 5 + 6! y- 7 + O (y~ 8 ) (19) 



4 Discussion 

In equations (fTUj) and (JT^J) the combination of terms r = aA£> e /<5 u0 is equal to l/z/ ( a_1 ) 
(from equation (JHJ) where f2„o = ^o)- Therefore r is the relaxation time for the charge on 
a dust grain with radius a" 1 to reach equilibrium with the ambient plasma, and V t r is the 
distance travelled by the test charge in this time. As remarked above, if we let a — > oo 
with p ~ a 3 the distribution h(a) tends to a monosized distribution with dust grain 
radius ao = (/9/3) 1//3 a _1 . Putting r = (/5/3) 1//3 r , the test charge response 4> c h given by 
equation (|13|) may be shown to reduce to the known results for a monosized distribution 
with a charge relaxation time tq 0111 ■ 
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